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We present preliminary results for the non-perturbative determination of the parameters in a 2+2 anisotropic 
action. The parameters are required to restore Lorentz invariance in long-distance correlators, using the static 
inter-quark potential. Comparison with analytical results is made and further applications are discussed. 



1. INTRODUCTION 

QCD on 3+1 anisotropic lattices |1I2I3| has 
proved to be a successful tool for the non- 
perturbative investigation of interesting physi- 
cal phenomema including the glueball spectrum, 
bb states, gluon strings, the spectral density. 
The method is successful for processes where the 
hadronic final state is static in the center of mass 
frame but it is perhaps less advantageous for pro- 
cesses in which the hadron momentum in the final 
state is non-negligible compared to its mass. In 
such cases, independent fine spatial directions for 
each momentum degree of freedom are needed to 
describe the system precisely. This increases dra- 
matically the computational cost and, when con- 
sidering an anisotropic approach, the difficulties 
linked with tuning the parameters. 

For processes where the final state hadron is 
strongly peaked at a momentum approximately 
equal to its mass, a 2+2 lattice discretization may 
provide a solution. It offers the possibility to ex- 
tend the range of final state momenta accessible 
to lattice calculations and enhance the overlap 
with experimental data. In addition, the study of 
high-momentum glueball form factors could ben- 
efit from such a program. 

A 2+2 discretisation of the gauge action was 
proposed in Refs | ll.~>| and the parameters in 
the action were determined at the one-loop level 
in perturbation theory. In this talk a non- 
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perturbative determination of these parameters 
is discussed. A 2+2 Wilson-like fermion action is 
discussed in Ref [7]. 

2. YANG-MILLS THEORY ON A 2+2 
LATTICE 

We start from the simple lxl plaquette action 

s g = PJ2 c ^[! - ^r Tr + (*))], (i) 

defined in a fixed hypercubic physical volume 
V = L 4 . The ultraviolet cutoff is larger in the 
z, t directions and £ = a c /a,f is the asymmetry 
ratio. In the following, the coarse and fine lattice 
spacings are denoted a c and a / respectively. The 
coefficients in Eq. ^ must be fixed by tuning 
their values so as to recover the correct continuum 
theory. This can be done either perturbatively [5J 
or non-perturbatively [5]. The perturbative ex- 
pansion is written 




e(l + V?)9 2 + 0(g*)) f-f; 

l + V$9 2 + 0(g*) c-f; (2) 

^(1 + r^V + 0( ff 4 )) c-c; 



and the calculation of the one-loop coefficients, 
rfpv was the subject of Ref [5]. We note in passing 
that the tree-level values of the are 

c cc = l/£ 2 , c cf = 1, and c ff = £ 2 . (3) 

A non-perturbative determination of these coeffi- 
cients is the focus of this talk. 
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We choose to fix the coefficients to restore ro- 
tational invariance for the static inter-quark po- 
tential. The procedure is, in principle, the same 
as for a 3+1 discretisation. However, for 3+1 any 
choice of the parameters leads to a well-defined, 
albeit unknown, scale ratio £. In this case, a non- 
perturbative determination of c M „ may be unnec- 
essary if, as is true for some applications, one 
is not interested in the physical mass-scale ratio 
between the coarse and fine directions. If such 
a determination is necessary, both the coarse- 
fine and coarse-coarse Wilson loops must be mea- 
sured. The same is true for the 2+2 case, where 
the potential between spatially separated color 
sources, measured using a fine temporal axis de- 
termines just one of the two free parameters in 
the action. The only drawback for the 2+2 dis- 
cretisation compared to the 3+1 is that such a 
determination must be done a priori. Rotational 
invariance of the static inter-quark potential must 
be imposed for all combinations of source separa- 
tion and transfer matrix axes. Alternatively, the 
slope of dispersion relations of physical (color sin- 
glet) states can be measured. This requirement 
can be observed in perturbation theory by exam- 
ining the static potential, given by 

V(R) = lim -i log < W(R,T) > . (4) 

Expanding Eq. ^ at order g A and handling the 
resulting Bessel functions for large T, it is easy 
to check that only the coefficients r)$ carrying 
one time index survive. A fine temporal axis will 
therefore give no information about r]^) , since 
constraining potential ratios in different direc- 
tions leads to a constraint only on c c //c// for 
2+2 and no constraint for 3+1. This constraint 
degeneracy is illustrated in Fig. Q] The solid and 
dashed lines represent the constraints on the spa- 
tial potential: firstly that V(n, 0, 0) = 7(0, 0, <fra) 
and secondly that V(n, n, 0) = V(n, 0, The 
outer solid and dotted lines are the one sigma er- 
rors from a fit to these constraints. Within the 
errors shown the lines are degenerate and cannot 
provide a non-perturbative determination of the 
c M „. Conversely, considering a coarse temporal 
axis results in no dependence on mj . For 2+2, 
only the ratio, c cc /c c / will be constrained. But 
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Figure 1. Degenerate constraint curves from the 
SU(3) spatial potential. 



now, this condition combined with the previous 
determination fixes the parameters uniquely. For 
3+1 this is the only meaningful condition. With 
a coarse temporal axis obtaining good fits can 
be more challenging. However, we find that with 
reasonable statistics good y^/d.o.f. is achieved. 
This non-perturbative determination of the 
is work in progress. 

3. PRELIMINARY RESULTS 

As a first step we show that the tuning de- 
scribed in the previous section is actually not as 
onerous as it might sound. Fig. [21 shows the po- 
tential in different directions, calculated using the 
tree-level coefficients. Even at the tree-level there 
is clearly no dramatic breaking of Lorentz invari- 
ance. This indicates that the tuning of param- 
eters is mild and the scheme can be of practical 
use. 

To demonstrate how powerful this scheme can 
be, Fig. shows the SU (3) potential calculated 
for £ = 6 on a 8 2 x 48 2 volume, in the range 
R = 0.04— 1.6/to. These results were obtained af- 
ter two days of workstation CPU time. The data 
show that the short distance term can already be 
calculated with very high precision. Longer runs 
and bigger volumes could easily allow one to reach 
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Figure 2. The S?7(3) static potential for £ = 6 
with tree-level coefficients c^ v , given in Eq. 01 
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Figure 3. The ST/ (3) static potential deter- 
mined with the one- loop values of the coefficients, 
c-ilv ■ The potential is measured over the range 
R = 0.04 — 1.6/m, for £ = 6. The errors on the 
points shown are smaller than the symbol size. 



even higher precision. The method is therefore 
immediately applicable to a very high resolution 
determination of the short-range inter-quark po- 
tential, which is of phenomcnological interest. 

4. OUTLOOK AND DEVELOPMENTS 

We have shown how 2+2 lattice QCD can be 
relevant to useful physics. Although parameter 
tuning is necessary, both the perturbative [5] and 
the non-perturbative determinations are feasible. 
The results for the static inter-quark potential 
demonstrate that this scheme can be a useful 
tool for exploring energy ranges higher than those 
currently available to lattice simulations. In the 
gauge sector the scheme can naturally be applied 
to a calculation of glueball form factors, and this 
is under investigation. An extension to the quark 
action allows further interesting phenomenologi- 
cal applications and is described in Ref [7]. 
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